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ABSTRACT 

We use numerical simulations of the time evolution of global textures to investigate the 
relationship between ordering dynamics and energy density in an expanding universe. 
Events in which individual textures become fully wound are rare. The energy density 
is dominated by the more numerous partially wound conhgurations, with median 
topological charge ctp ~ 0.44. This verihes the recent supposition [If] that such 
partially wound conhgurations should dominate the cosmic microwave background. 
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I. Introduction 


Recent observations [1, 2] of anisotropy in the cosmic microwave background (CMB) 
provide strong evidence for primordial density perturbations. Several models predict 
the statistical distribution of these perturbations. In inflationary models, density 
perturbations result from parametric amplihcation of quantum fluctuations during 
the inflationary epoch [3], while in topological defect models the density perturbations 
result from large-scale ordering effects following a phase transition with broken 
symmetry [4, 5, 6, 7]. 

Observational tests of these models are complicated by the fact that 
cosmological models predict the CMB anisotropy to be a single realization of a 
stochastic process whose properties (power spectrum, phase correlations, etc.) are 
predicted only for an ensemble average of equivalent realizations. Any single 
realization may be expected to vary about the ensemble average. Although a large 
body of work exists providing analytic derivation of ensemble average properties of 
inflationary models, the non-linear nature of defect models preclude exact analytic 
results. Tests of defect models must use either time-consuming numerical simulations 
or an approximate analytic model. 

Global textures are one class of defect models for which an approximate 
analytic model has been proposed [7]; however, the predictions of this approximate 
model do not agree with the results of a complete simulation [8]. One of the main 
assumptions of this model is that the CMB anisotropy is dominated by the most 
energetic isolated defects. We use a numerical simulation of the complete theory to 
investigate this assumption and hud that the dominant contribution comes instead 
from the far more numerous but less energetic defects. This was hrst suggested in [If] 
for a non-expanding universe. By doing the simulation for an expanding universe, we 
demonstrate that the less energetic defects must be considered in any model of the 
texture generated CMB anisotropy. 

Topological defects such as monopoles, domain walls, or cosmic strings involve 
regions of space where the order held remains in the higher-energy symmetric state. 
The energy density in these models is dominated by the order held potential. Textures 
[6] are a class of topological defects with the property that the broken symmetric phase 
vacuum manifold of the order held is of the same dimension as the spatial geometry. 
This allows the order held always to stay on the vacuum manifold, regardless of initial 
conditions. Perturbations in the energy density are driven solely by the order held’s 
kinetic energy. 

The topological charge of a region of space is dehned as the fraction of the 
vacuum manifold covered as the order held is mapped into the manifold [9]. At the 
time of the phase transition, the order held in each causal domain will randomly 
fall into any of the degenerate vacua. As the different domains come into causal 
contact, the held may develop a non-trivial topological charge, forming a wound 
conhguration or “knot”. Derrick’s theorem [10] states that the spatial extent of any 
region of non-zero gradient energy must shrink, causing the order held to become 
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increasingly tightly wound around the vacuum manifold. The kinetic energy in this 
wound conhguration will eventually become greater than the potential barrier to the 
symmetric state, lifting the order held from the vacuum manifold and allowing the 
texture to “unwind” as it reduces the topological charge. 

Since this unwinding process is so energetic, it is often viewed as the dominant 
contribution to the energy density. Turok & Spergel [7] provide an analytic model of 
the CMB anisotropy resulting from a distribution of unwinding spherical knots. The 
results of this model are not reproduced by numerical simulations [8]. On the other 
hand, Borrill et. al. have pointed out that unwinding events are rare: they found 
only 40 such events in 1000 simulations [9]. This agrees with analytic predictions 
[6] (see also [4, 15]). They also hud partially-wound texture conhgurations can still 
generate anisotropies similar to those generated by events which do unwind [11]. 

We have simulated the order held dynamics in the non-linear a model 
approximation to investigate the distribution of gradient energy density. We identify 
both the wound and partially-wound events in our simulation. The unwinding events 
are indeed rare: we hud only 3 unwindings events out of over 400 texture events. 
Although unwindings are the most energetic individual events, most of the energy 
density is associated with the more numerous partially-wound texture conhgurations. 


II. Analysis 

Our calculations are done in a spatial hat Robertson-Walker homogeneous universe, 
with metric given by 

ds — —dt T (^dx T dy T dz ) = (2(t) (— dx T dx T dy T dz ). (1) 

We consider the evolution of a real global held a = 1,. . . , 4, whose dynamics are 
specihed by the Lagrangian density 

^ ( 2 ) 

The function is such that the held $ undergoes a phase transition at T = Tc 

from a symmetric phase (T > Tc) to a broken symmetry phase (T < Tc). The 
symmetric phase corresponds to 4>‘^{T) < 0, giving the held an 0(4) symmetry. When 
<t>‘^{T) > 0, the held falls into the broken symmetry phase with 0(3) symmetry. The 
vacuum manifold becomes = 0(4)/0(3). For the range of dynamics we consider, 
we assume T Tc and <f>‘^{T) = do? ^ constant. 

Textures form as horizons meet, which has the length scale of Hubble’s constant 
H. The natural length scale for textures is the inverse of the radial excitation mass of 
the held. It is at this scale that we expect a collapsing “knotted” texture to unwind 
[12]. Borrill [13] has pointed out that for a GUT scale phase transition, these scales 
differ by 10®° at the time of decoupling. Any attempt to numerically simulate this 
model can not hope to cover this range of scales and an approximation must be made. 
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In this work, we follow [8] and use the non-linear a model (NLSM). This takes the 
held’s mass to inhnity, forcing $ to always be on the vacuum manifold. 

The NLSM model is obtained when the Lagrangian density is rewritten as 

= (3) 

where the constraint that the held stay on the vacuum manifold is now imposed via 
the Lagrange multiplier A (the held rescaling —)• $a/^o is assumed). Variation of 
this Lagrangian leads to the equations of motion 


= [|V$|2 - |$|2] (4) 

To investigate this theory via computer simulation, we simulate the order held 
evolution on a discrete cubic lattice. We map $“(x,f) —)• where i labels the 

spatial cells and n is the temporal index. Without loss of generality, we take the grid 
spacing 6x = 6y = 6y = 1 and set the scale via our choice of the time interval St. 
Varying a discrete version of Eq. (4) [8], we get the equations of motion 


$ 


a 




(5) 


where 

denotes a^(T„ ± is chosen to enforce the constraint 

Squaring Eq. (5) and solving for A gives: 


( 6 ) 

1 . 


A, = ± (7) 

Expect for the case mentioned below, we always take the positive square root. 

The NLSM implementation provides no method by which a knotted 
conhguration unwinds. Since $ is constrained to the vacuum manifold, the equations 
of motion no longer allow for unwinding. In effect, we have assumed that the spatial 
scale for unwindings is smaller than the inter-cell distance. If we can recognize a 
wound conhguration once it has collapsed to be comparable to the cell size, we may 
add an unwinding mechanism “by hand”. Eor each spatial cell, we dehne an alignment 

= ( 8 ) 


where the sum on j runs over the cell’s six nearest neighbors. A cell is counter- 
aligned if ai < 0, when it has become completely wound on the scale of the cell 
spacing. This can be thought of as a situation where the topological charge within a 
3x3x3 cube centered around the counter-aligned cell becomes > 0.5. The unwinding 
then proceeds by choosing the negative sign in Eq. (7): the cell is moved across the 
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vacuum manifold to a configuration that removes the topological charge [8]. Eq. (7) 
is thus modihed to become 


A* = + sign(a*)^l - (9) 

and the knotted textures are unwound once they shrink below the cell grid scale. The 
energy density at each point is given by the rr-component of the stress energy tensor: 

P = \\^? + ( 10 ) 

and we use the obvious discretised version in our analysis. 


III. Numerical Results 

We use the above algorithm to investigate the distribution of the energy density during 
the held evolution from the time of decoupling onward. Since the decoupling time is 
long after the time of the phase transition, we consider held evolution starting from 
an initially smooth held conhguration. This is achieved by initializing each point on 
the simulation grid to be uniformly distributed on the vacuum manifold and allowing 
the equations of motion to evolve the held until it “settles down”. Initially, the ads 
will be uniformly distributed in [—1,1]. If the unwinding mechanism (Eq. (9) instead 
of Eq. (7)) were used, it would lead to an unstable held evolution. Instead, we use Eq. 
(7), choosing strictly the positive square root and thus keeping the gradients small. 
We continue until the total number of counter-aligned cells stops decreasing. Then 
the unwinding mechanism is enabled and Eq. (9) is used to evolve the held. This 
procedure allows the held evolution algorithm to generate the required smooth initial 
conditions. The late time behaviour of the held is insensitive to the precise time that 
the unwinding is enabled. 

We have run a simulation of a 64^ element spatial cube for 32 conformal time 
units (c= 1, 8x = l). Since periodic boundary conditions are used, this is the maximum 
length of time the simulation can run for and still be viewed as physical. We took 
8t = 0.2. The evolution was insensitive to taking 8t —)• 8t l2. The evolution was also 
insensitive to using a 32^ elements spatial cube. The simulation became stable for 
T > 10, at which time the unwinding was enabled. 

The simulation shows the expected behaviour as regions of above average 
energy density form, shrink in spatial extent while the energy density rises, and then 
re-expand while the energy density falls back to the average background level. The 
time scale for this rise and fall is of order At ~ 3. We call such a process an event. 
Typically we hnd the energy density at the center of an event to be at least 4(7 above 
the average. 

Eig. 1 shows the energy density vs cell alignment for the 20 most energetic 
events. The alignment decreases and the energy increases as the texture becomes 
increasingly wound. There were three unwinding events [a < 0) in this simulation: 
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the energy density increases by an order of magnitude as the knotted texture unwinds. 
The rest of the events never become completely wound. The alignment decreases to a 
minimum a > 0 and then starts to increase again. Shortly after the alignment stops 
decreasing the energy density stops increasing and the event dissipates. 

Fig. 2a shows the distribution of minimum alignments for all events. The 
distribution’s median is a ~ 0.48. The most likely event is one where the texture 
never becomes completely knotted. Although the energy density in these partially- 
wound events is an order of magnitude below the fully knotted conhguration, they 
are so numerous as to dominate the total energy density of the simulation. Fig. 2b 
shows the energy density integrated over time for all events in the simulation vs the 
minimum alignment of each event. This distribution’s median is ctp ~ 0.44. The 
dominant energy contribution in a texture cosmology comes from partially wound 
conhgurations, with only minor contributions from unwinding knots. 


IV. Discussion 

The Sachs-Wolfe effect [14] relates CMB anisotropy to perturbations in the 
gravitational potential, and is usually given in terms of an integral over the the time 
rate of change of the gravitational potential. For texture models, the gravitational 
potential results from the kinetic energy density of the order held. There does not 
appear to be a great variation in the spatial prohle of the individual events; hence, 
the contribution of each texture event to the CMB anisotropy should be proportional 
to the integrated energy density in the event. 

Numerical simulations of the time evolution of global textures demonstrate 
that the integrated energy density is dominated by partially-wound conhgurations 
and not by unwinding knots. The gradient energy density due to texture events is 
centered about textures with alignment tip ~ 0.44. This result proves Borrill et. al. 
[If] claim, but now in the more realistic scenario of an expanding universe. 

Ffforts to obtain observable predictions for specihc texture models have 
either relied on computationally expensive numerical simulations or on simplifying 
assumptions. The hrst method suffers from logistical problems in that proper 
assessment of cosmic variance requires hundreds to thousands of realizations for each 
choice of model parameters. The second method typically suffers from assumptions 
(spherical symmetry, uncorrelated placement of subsequent textures, fully wound 
textures) that may not be well founded. We have demonstrated that any simplihed 
texture model must consider the full distribution of partially-wound events along 
with the rarer high-energy events. We are currently attempting to derive a 
statistical description (distribution of shapes, spatial and temporal correlations) 
of the gravitational potential in texture cosmologies. If a sufficiently compact 
description can be obtained which retains the distinctive features of a full-blown 
numerical simulation, signihcant computational savings may be realized in Monte 
Carlo simulation of CMB anisotropy for texture models. These techniques may 
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then be incorporated into other analytic models [15, 16] to compare specihc texture 
cosmologies with the observed CMB anisotropy on various angular scales. 


Acknowledgements 

We would like to thank Shoba Veeraraghavan for her time and helpful discussions. 
One of us (NGP) would also like to thank Fred Lang of Catholic University of America 
for making the funds and resources available to facilitate this work. 


7 



References 

[1] G. F. Smoot et al.^ Astrophys. J. 396 LI (1992). 

[2] C. L. Bennett et al.^ Astrophys. J. 436 423 (1994). 

[3] J. M. Bardeen, P. J. Steinhardt and M. S. Turner, Phys. Rev. D, 28 679 (1983); 
A. Guth and Y-S. Pi, Phys. Rev. Lett. 49 1110 (1982); S. Hawking, Phys. Lett. 
115B 295 (1982); A. A. Starobinskii, Phys. Lett. 117B 175 (1982). 

[4] A. Vilenkin and E. P. S. Shellard, Cosmic Strings and Other Topological Defects 
(Cambridge University Press, Cambridge, 1995). 

[5] D. P. Bennett and S. Rhie, Astrophys. J. 414 L7 (1993) (global monopoles); A. 
Vilenkin, Phys. Rep. 121 263 (1965); A. Stebbins et. a/., Astrophys. J. 322 1 (1987) 
(cosmic strings); A. Stebbins and M. S. Turner, Astrophys. J. 329 L13 (1989); M. 
S. Turner, R. Watkins and L. M. Window, Astrophys. J. 367 L43 (1991) (domain 
walls). 

[6] N. Turok, Phys. Rev. Lett. 63, 2625 (1989). 

[7] N. Turok and D. Spergel, Phys. Rev. Lett. 64 2736, (1990). 

[8] U.-L. Pen, D. Spergel, and N. Turok, Phys. Rev. D 49 692 (1994). 

[9] J. Borrill, E. J. Copeland, and A. R. Liddle, Phys. Rev. D 47, 4292 (1993). 

[10] G. H. Derrick, J. Math. Phys., 5 1252 (1964). 

[11] J. Borrill, et. al., Phys. Rev. D 50 2469 (1994). 

[12] N. Turok, Phys. Scr. T36, 135 (1991) 

[13] J. Borrill, Phys. Rev. D 50, 3676 (1994). 

[14] R. Sachs and A. Wolfe, Astrophys. J. 147, 73 (1967). 

[15] J. C. R. Magueijo, Analytical Models of large-angle CMBR anisotropies from 
textures, MRAO/1792, DAMTP/94-101, to be published in Phys. Rev. D. 

[16] R. Durrer, A. Howard, and Z.H. Zhou, Phys. Rev. D 49, 681 (1994); R. Durrer 
and Z.H. Zhou, Phys. Rev. Lett. 74 1701 (1995). 


8 




Alignment 


Figure 1: Evolution of texture events viewed in the energy density-alignment plan. 
Diamonds denote points in the evolution where the alignment was negative (an 
unwinding event). 
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Figure 2: (a) Distribution of minimum alignments of event histories, (b) Energy 
density (arbitrary units) integrated over time, sorted by the minimum alignment of 
each event. The median for energy distribution is tip ~ 0.44, from texture events that 
are only partially wound. 
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